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In approximate dynamical equations, inhomogenous classical (mean) gauge and Higgs fields are coupled to 
quantized fermions. The equations are solved numerically on a spacetime lattice. The fermions appear to equili- 
brate according to the Fermi-Dirac distribution with time-dependent temperature and chemical potential. 



1. The real-time path integral for quantum 
fields is very difficult to evaluate numerically and 
approximations need to be made before giving 
the problem to the computer. Two types of ap- 
proximations are currently in use: classical and 
gaussian, such as Hartree, large N. The clas- 
sical approximation gives valuable nonperturba- 
tive results M but suffers complications due to 
(Rayleigh-Jeans type) divergencies. The gaussian 
approximation has the benefit of staying within 
the quantum domain where we know how to deal 
with divergencies, but it is not good enough for 
large times. Naturally, one would like to com- 
bine the good aspects of both approximations 
(21. A crucial test is to see whether the sys- 
tem equilibrates quantum-like, and not classical 
equipartition-like, despite the fact that one is just 
solving a large number of coupled nonlinear equa- 
tions which conserve energy. This is one motiva- 
tion for the present study. Another is the intrinsic 
interest in the complicated nonperturbative dy- 
namics of the abelian Higgs model with fermions. 

2. The 1+1 D abelian Higgs model coupled 
axially to fermions is qualitatively similar to the 
electroweak sector of the Standard Model. As for 
the SU(2) case in 4D, it can be rewritten in a form 
with vectorial gauge couplings and Majorana- 
Yukawa couplings. For N — > oo fermion replicas, 
the equations of motion reduce to a classical field 
approximation for the bosonic variables, with a 



quantal fermion backreaction ||: 

d„F^ v + e 2 i{D"ip* ip - v*D u tp) 

+ (e72)(^ 7 » = 



0. 

0. 



where we have specialized to zero Yukawa cou- 
pling. The fermion backreaction is specified as 
follows. Introduce a complete set of orthonormal 
mode functions u a , v a for the fermions, which sat- 
isfy the Dirac equation 

Next define the fermion operator ip, 
i>(x) = ^2[b a u a (x) + d) a v a {x)], 
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in terms of annihilation and creation opera- 
tors b a , &T, , .... The fermion back reaction is 
then specified by the initial conditions (b\b a i) = 

$olol' > (^a^a') &aa' ) etc. 

The above system of equations has been im- 
plemented on a lattice using Wilson's fermion 
method for the spatial derivative and the stag- 
gered fermion interpretation for a 'naive' dis- 
crete time derivative [||. Usual expectations 
on fermion number non-conservation tied to 
sphaleron transitions are correctly represented on 
the lattice. For simplicity we continue with con- 
tinuum notation. 

3. Here we are especially interested in ther- 
malization properties of the fermions. We tested 
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Figure 1. Chcrn-Simons number C and axial 
charge (Qs)/^/ versus time in units 1/e. 



Figure 2. Least squares straight line fits to 
ln(N~ vp — 1) versus ap. 



|I| for local (in time) equilibration of fermions 
by comparing their distribution function with the 
Fermi-Dirac distribution. The distribution func- 
tion was identified from the equal-time fermion 
two point function, averaged over space (a circle 
with circumference L), 

1 f L 

S(z,t) — — J dx \ip(x, t)ip(x + z, t))g f . 

Here g.f. indicates a complete gauge fixing. Alter- 
natively, the two point function can be rendered 
gauge invariant by supplying a parallel trans- 
porter U(x,y) — exp[ — i J y dz Ax{z)/2). We have 
used the latter method, but it is actually closely 
related to complete Coulomb gauge fixing B. If 
the fermions were free, the Fourier transform 



>(*) = y^W + 



dze~ lpz S(z,t) 



would be given in terms of distribution functions 
N p , N p als follows: 



Tr Sfat) = [1 - N p (t) - JV_„(t)] 
TttfSfat) - [l-N p (t)-N- p {t)} 
Tri 7 °S(p,t) - 1 — N p (t) + N- P (t), 



m p (t) 

LOp(t) ' 

P 

U! p (t) ' 



with Tr75S'(p, t) = because of parity invari- 
ance. For free fermions N p , N p and m are 
time-independent. Assuming that the interacting 
model can be described approximately by quasi- 
particles, we now use the above equations to de- 
fine N p (t), N p (t) and m p (t). In a non-equilibrium 
situation they will depend on time. 

The simulations had the following parameters: 
rif = 2 flavors (related to fermion doubling in 
time), spatial size m v L « 6.4, A/e 2 = 0.25 
{ttiaL w 9), coupling e 2 /m 2 w 0.25, with spa- 
tial lattice spacing am v w 0.10, temporal spacing 
ao/a = 0.005 and L/a = 64 spatial lattice sites. 

Fig. |l| shows the Chern-Simons number C = 
— J dxAi/2-K and the axial charge (Qs)/n/ for a 
simulation starting with a fcrmionic vacuum (i.e. 
n a ,n a = 0) and some kinetic energy stored in 
a few low momentum modes of the Higgs field. 
The anomalous fermion number non-conservation 
equation A(Q$) = n/AC, is well obeyed since the 
two curves are indistinguishable (initially (Q5) = 
C = Q). The oscillations correspond roughly 
to the basic period 2ir/rriA- To smoothen these 
we average S(p, t) over a time interval t av before 
extracting the distribution functions. We used 
ei av = 4 and studied the behavior of iV avp (i) = 
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Figure 3. Effective temperatures T(t) — l//3(£) 
versus time obtained from fits as in figure 2. 



Figure 4. As in figure 3 for the effective chemical 
potential. 



[N p (t) + N- p (t)]/2. As expected, m p (i) is effec- 
tively zero. Surprisingly, N avp (t) resembles quite 
fast a Fermi-Dirac distribution for inverse tem- 
perature (3 and Qs-chemical potential [i: 



{Q 5 )/L 



iffx/ir. 



E n — 



{exp[f3(E p - /j,q 5p )} 
\p\, Q5 P =p/\p\ 
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(the axial charge of a fermion depends on the 
sign of p). Fig. 2 shows ln(A^p — 1) versus 
ap at various times. We see linear behavior, 
M N avp - 1) « P{t)[\p\ ± suggesting local 

(in time) equilibrium. The distribution functions 
of modes with ap > 0.5 are consistent with zero, 
so these modes are practically not excited. Fur- 
thermore, since the relevant p are small in lattice 
units, discretization effects are reasonably small. 
To achieve this the initial energy stored in the 
Bose fields has to be in the relatively low momen- 
tum modes only - the fields are far from classical 
equilibrium. 

Figs. H and 4 show the effective temperature 
and chemical potential as a function of time. Note 
that T p> o(t) ~ T p< o(t). If the fermions were free, 
then their energy and axial charge densities would 
follow from the Fermi-Dirac distribution accord- 
ing to 



Conversely, Ep and (Q5) imply an effective tem- 
perature and chemical potential; these are also 
plotted in Figs. ||[4 (data labeled Ep,C resp. 
C). The Ep, C-temperature appears systemati- 
cally lower that that from A avp . This may be 
due to the fact that the fermions are not free. 

4. In conclusion, we see evidence for fast equi- 
libration of fermions coupled to classical Bose 
fields. It is important that the classical fields 
can be spatially inhomogeous, since this allows 
the fermions to scatter nontrivially, e.g. by their 
(screened) Coulomb interaction. The Bose fields 
have created fermions and lost some energy, and 
they are not in (classical) equilibrium in the time 
span shown here. More details on this work can 
be found in Q. 
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